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ABSTRACT 

Context. Recent results by Bensby and collaborators on the ages of microlensed dwarf and subgiant stars in the Galactic bulge have 
challenged the picture of an exclusively old stellar population, because ages significantly younger than 9 Gyr have been found. 

Aims. However, these age estimates have not been independently confirmed with different techniques and theoretical stellar models. 
One of the aims of this paper is to verify these results by means of a grid-based method. We also quantify the systematic biases 
that might be induced by some assumptions adopted to compute stellar models. In particular, we explore the impact of increasing 
the initial helium abundance, neglecting the element microscopic diffusion, and changing the mixing-length calibration in theoretical 
stellar track computations. 

Methods. We adopt the SCEPtER pipeline with a newly computed stellar model grid for metallicities [Fe/H] from -2.00 dex to 0.55 
dex, and masses in the range [0.60; 1.60] Af 0 from the ZAMS to the helium flash at the red giant branch tip. By means of Monte 
Carlo simulations we show for the considered evolutionary phases that our technique provides unbiased age estimates. 

Results. Our age results are in good agreement with Bensby and collaborators findings and show 16 stars younger than 5 Gyr and 
28 younger than 9 Gyr over a sample of 58. The effect of a helium enhancement as large as AK/AZ = 5 is quite modest, resulting in 
a mean age increase of metal rich stars of 0.6 Gyr. Even simultaneously adopting a high helium content and the upper values of age 
estimates, there is evidence of 4 stars younger than 5 Gyr and 15 younger than 9 Gyr. For stars younger than 5 Gyr, the use of stellar 
models computed by neglecting microscopic diffusion or by assuming a super-solar mixing-length value leads to a mean increase in 
the age estimates of about 0.4 Gyr and 0.5 Gyr respectively. Even considering the upper values for the age estimates, there are four 
stars estimated younger than 5 Gyr is in both cases. Thus, the assessment of a sizeable fraction of young stars among the microlensed 
sample in the Galactic bulge appears robust. 
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1. Introduction 


Recently a debate on the age of the stars in the Galactic bulge has 
developed. Until few years ago, the picture of a Galactic bulge 
composed of old stars was commonly accepted. The evidence 
from Hubble Space Telescope photometry of the Galactic bulge 
main-sequence turnoff points toward an exclusively old stellar 
population, with ages greater than about 10 Gvr (e.g. Zoccali 

2003; Brown et al.l2010,: Clarkson et alj201 1; Valenti et all 


et al. 


201311 . The similarity between the Galactic bulge and Galactic 
globular cIu sters suggested an age estimate of 12-13 Gvr ( Marin- 
Franch et al. l2009 l) 

However, in recent years a different scenario coming from 
gravitational microlensing results was proposed. iBensbv et alJ 
( 201 1. 2013! ) obtained high-resolution spectra of 58 bulge main- 
sequence and subgiant stars, from which they derived measure¬ 
ments of T e fj, log#, and [Fe/H]. Surprisingly, the comparison of 
the observations with theoretical isochrones revealed a signifi¬ 
cant fraction of young and intermediate-age stars. In fact, the 
IBensbv et al.l J2013h results showed the presence of 13 stars out 
of 58 (22%) younger than 5 Gyr and 31 stars (53%) younger than 
9 Gyr. 

As a consequence of these results, there has been some ef¬ 
fort in the literature to reconcile the two pictures. A first attempt 
was made bv lNataf & Gouldl ( )2012 ). who suggested that taking 
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into account the possible Galactic bulge helium enhancement, 
the differences between the two sets of age estimates should van¬ 
ish. However, this result was questioned bv lBensbv et al.l d2013l) . 
who claimed that the effect of the helium enhancement is insuf¬ 
ficient to contradict the conclusion of a substantial presence of 
stars at young and intermediate ages. Furthermore, a possible 
contamination of the microlensing s amp le by disk star and bias 
toward a lower age was proposed (lNatafll20~13h . However, none 
of these distortions seem to be large enough to refute the Bensby 
et al. (120131) ' findings. More recently. iNess et al.l (120141) sug¬ 
gested that the discrepancy could be due to the different spatial 
distributions in the Galactic bulge of old and young stars. 

Given the importance of the age estimates by IBensbv et al.l 
(1201 3h i n the framework o f Galactic bulge formation scenarios 
(see e.g. INess et aill2()l4l and references therein), we feel that 
an independent analysis of these results - adopting a different 
estimation technique and different theoretical stellar models - 
and an analysis of potential systematic biases is mandatory. The 
purpose of the present work is to provide such an independent 
check adopting the SCEPtER pipeline (IValle et alJl2014bllah . 

We also explore in detail the impact on these age estimates of 
changing some parameters or prescriptions required to compute 
a grid of stellar models that are not yet tightly constrained. In 
particular, we analysed the effect of a helium-rich Galactic bulge 
by computing a grid of stellar models with a helium-to-metal 
enrichment ratio AT/AZ much higher than the standard one (i.e. 
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AY/AZ = 5 instead of 2). Moreover, we checked the systematic 
differences in age estimates due to the inclusion or not of the 
microscopic element diffusion in stellar evolution calculations, 
whic h is an important bias sou rce in gri d-based techniques ( Va11e 
et al. l2014bHal) . Finally, we quantified the effect of changing the 
mixing-length calibration in stellar models. 

The structure of the paper is the following. In Sect. [2] we 
discuss the estimation method and the stellar-model grids used 
in the process. The age estimates are presented in Sect. [3] and 
the analysis of the impact of changing the grid of models by tak¬ 
ing into account a much higher helium abundance, by neglecting 
the element microscopic diffusion, and by increasing the mixing- 
length parameter is performed in Sect. [4] Some concluding re¬ 
marks can be found in Sect. [5] Finally, a Monte Carlo validation 
of the adopted method for stars in evolutionary stages similar to 
those of the microlensed sample is presented in Appendix [A] 


2. Methods 


The age estimates obtained in this paper are based on our re¬ 
cently developed pipeline SCEPtER, which is extensively de¬ 
scribed in (IValle et al.l2014bl lah and briefly summarized here. We 
let ^ be a star for which the following vector of observed quanti¬ 
ties is available: q s = {7' e fr.^, [Fe/H]^, logg^}. Then we let cr = 
{o"(7’ e ff 5 ),cr([Fe/H] 5 ),cr(logg 5 )} be the nominal uncertainty in 
the observed quantities. For each point j on the estimation grid 
of stellar models, we define q■' = {T e sj, [Fe/H| ; , logg ; ). We let 
J/ j be the likelihood function defined as 


A' = 

where 


fl^— 

' = | V27TO-/ 


x exp I — — 1, 


X 


-2 

i= 1 


qr-Vi 


v cr ' , 


(i) 


( 2 ) 


The likelihood function is evaluated for each grid point 
within 3cr of all the variables from S', we let .£ max be the maxi¬ 
mum value obtained in this step. The estimated stellar mass, ra¬ 
dius, and age are obtained by averaging the corresponding quan¬ 
tity of all the models with likelihood greater than 0.95 x _£ max . 

The technique allows a Monte Carlo confidence interval to be 
constructed for the estimates. To this purpose a synthetic sample 
of n stars is generated, following a multivariate normal distribu¬ 
tion with vector of mean q s and covariance matrix X = diag(cr). 
A value of n = 10000 is usually adopted since it provides a 
fair balance between computation time and the accuracy of the 
results. The medians of the estimated mass and age for th n 
objects are taken as the best estimate of the true values; the 16th 
and 84th quantiles of the n values are adopted as a 1 cr confidence 
interval. 


2.1. Standard stellar models grid 


The standard estimation grid of stellar models was obtained 
using the FRANEC stellar evolution code (Degl’Innocenti 


et al. 2008 ; Tognelli et alJboilh , as it was adopted to com- 
pute the Pisa Stellar Evolut i on Data BasJ] for low-mass s tars 

(jPeH’Ornodarme et al.ll2012i : IPelTOmodarme & Vallell201 3). 

The grid consists of stellar models in evolutionary stages 
from the zero-age main sequence (ZAMS) to the helium flash 
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at the tip of the red giant branch (RGB). Models of age greater 
than 15 Gyr were excluded from the grid. Models were com¬ 
puted for stellar masses in the range [0.60; 1.60] M e with a step 
of 0.02 Mr.,. The initial metallicity [Fe/H] was assumed in the 
range [-2.0; 0.55] with a step of 0.10 dex from -2.0 to -1.0 
and 0.05 dex for [Fe/H] > -1.0. The solar scaled heavy-element 
mixture bv lAsnlund et al . 6200 90 was adopted. For [Fe/H] lower 
than -0.5 dex we accounted for an a-elements enhancement 
[ar/Fe] = +0.3. The initial helium abundance was obtained us¬ 
ing the linear relation Y — Y p + ~^Z. The primordial 4 He abun¬ 
dance valueT a _^0.2485 from WMAP was adopted (Cyburt 


et al. l2004t ' SteiemaiJ 20061; iPeimbert et al.l |2007al b|)~ a ssum¬ 
ing A Y/AZ = 2 (Pagel & Portinarif 19981 : 1.Timenez et al.ll2003l : 
iGennaro et al.ll2Q10l) . The models were computed assuming our 
solar-scaled mixing-length parameter cr,,,! = 1.74. Atomic diffu¬ 
sion was included adopting the coefficients given bv lThoul et al.l 
(119941) for gravitational settling and thermal diffusion. Convec¬ 
tive core overshooting was not taken into account. Further de¬ 
tails on the input of the code can be found in 1 Valle et al j 6 2009 1 
1201 3 alibi) . 

A Monte Carlo estimation of the expected errors and pos¬ 
sible biases of our technique for stars similar to those in the 
microlensed bulge sample, adopting the available set of obser¬ 
vational constraints, is reported in Appendix lAl The analysis 
showed that the method is unbiased for the considered evolu¬ 
tionary phases. 

The adopted estimation technique and the grid of stellar 
mod els differ under many aspects from those used by Bensbv 
et al. (120131) , therefore a confirmation of their finding would 
add conside rable robustness to the age results. iBensbv et alJ 
(120111120131) adopted a maximum likelihood functional maxi¬ 
mization, where the errors on the estimates are derived from the 
likelihood p rofiles. They used the stellar models by Yonsei-Yale 
isochrones (IDemaraue et al.ll2004 ). which rely on input physics 
different from our own reference. In particular, they adopt the 
solar mixture from Grevesse et al. ( 199 6). the low temperatures 


opacities by Alexander & Ferguson (1994), the conductive opac¬ 
ities by Hubb ard & Lam pel (1 19 691) , the O PAL equation of state 
in the 1996 version (Rogers et all 1 19961) , helium but no metal 
micro scopic diffusion, and the cross sections listed in Bahcall 
(119891) , Moreover, they take into account a mild convective core 
overshooting following the recipe described in IDemaraue et all 
(120041) . 


3. Age estimates 

The observational data - with their errors - of the 58 microlensed 
bulge stars jitudied in this paper have been obtained from Table 4 
in IBensbv et al.l (120131) . For each object, we estimated the mass 
and the age by means of the SCEPtER technique, as detailed 
in Sect. [2] Figure |T| and Tab. Q] show our results together with 
the observed [Fe/H]. Figure [2] shows the comparison between 
our age estimates and those bv iBensbv et al . 12013?) (left panel), 
and the age differences between the former and the latter as a 
function of the metallicity (right panel). The error bars on the 
age differences are obtained by adding in quadrature the errors 
provided by the two techniques. As one can see, the two inde¬ 
pendent results generally agree with each other, as it is apparent 
that all the differences are consistent with zero. A comparison 
of the two sets by means of a paired r-test showed a smafl mean 
difference of -0.08 Gyr (95% confidence interval [-0.46; 0.30] 
Gyr), visible in Fig. [2] 

This agreement is further confirmed by comparing the em¬ 
pirical cumulative distributions of stellar ages estimated by the 
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Fig. 1 . [Fe/H] versus SCEPtER age estimates for the 58 microlensed 
stars. 



age (Gyr) 


Fig. 3. Empirical cumulative distribution functions for the 58 es¬ 
timated ages obtained by SCEPtER (black line) and by Bensbv et alJ 
d2013h (red line). 


two techniques, shown in Fig. [3] A two-sample Kolmogorov- 
Smirnov test did not reveal significant differences between the 
two empirical distributions (Z) = 0.105, p-value = 0.92) im¬ 
plying that the age distribution can be considered consistent be¬ 
tween the two estimation methods. 

It is apparent from Fig. [3] that both the techniques indepen¬ 
dently estimate a sizeable fraction of young stars. We found 16 
stars younger than 5 Gyr and 28 younger than 9 Gyr that is, 28% 
(95% confidence interval [17%; 41%]) and 48% (95% confi¬ 
dence interval [35%; 62%]) of the sample, respectively, whereas 
iBensbv et al. (2013) found 13 (22%) and 31 (58%). 


4. Effect of changing the grid of stellar models 


Basu et al 


The determination of stellar ages necessarily requires theoretical 
stellar models, which in turn depend on prescriptions and pa¬ 
rameters that have not yet been tightly constrained (tValle et al.l 
1201335) . As a consequence, the resulting age estimates are af¬ 
fected by non-negligible uncertainties and, possibly, biases. In 
recent years, several studies have been devoted to the analysis of 
the performance of grid based techniq ues (e.g. iGai et al.ll2011l : 

[2012; Mathur et al .1120 121) . In this perspective, in 
T2014hlU we carried out a detailed statistical anal¬ 
ysis of the main biases affecting any grid-based recovery pro¬ 
cedure. In this section we will discuss the impact on the age 
estimate of the 58 microlensed bulge stars of increasing the ini¬ 
tial helium abundance, switching off microscopic diffusion, and 
changing the mixing-length calibration adopted to compute the 
stellar models. The main aim of the test is to check the robust¬ 
ness of the assessment of a sizeable fraction of young bulge stars. 


Valle et al 


4.1. Effect of increasing the initial helium abundance 


Nataf & Gould (20 


2) suggested that the discrepancy in age be- 


_^lisuggi 

tween IBensbv et al] (120131) results and those inferred by photo¬ 
metric data would vanish whenever an enhancement in the he¬ 
lium content of the bulge stars was addressed. iNataf & Gouldl 
(120121) . studying the effect on theoretical isochrones of helium 
abundance variations, discussed the systematic effects on age es¬ 
timates that could arise when observational data were analysed 
with standard instead of helium-enhanced isochrones. In par¬ 
ticular, they found that the predicted trend of the difference be¬ 
tween the true absolute magnitude and the spectroscopically in¬ 
ferred one should be positively correlated with the spectroscopi¬ 
cally inferred stellar mass, a trend that they found also present in 
the IBensbv et al. ( 201 1) data. Thus they claimed that a value of 
helium-to-metal enrichment ratio of AT/AZ » 5 would reconcile 
the two sets of age estimates. In contrast, IBensbv et al] (12013l) 
showed that the discussed trend is no longer present in the ex- 
tcnde d stellar sample; they rejecte d the hy pothesis bv Nataf & 
Gould (i2012l ) and concluded that a sizeable sample of young and 
intermediate age stars would still be present even assu ming high 
helium abund a nces. However, neither Nata f & Gouldl (1 201 2l) nor 
IBensbv et al.l (1201 3l) performed a statistically detailed compar¬ 
ison of the age estimates obtained with different values of the 
original helium abundance. Moreover, the impact of assuming 
a different initial helium value is strongly dependent on the ob¬ 
servational quantities used to constrain the ages. As an example, 
in I Valle et al.l (l2014al) we have recently shown that in the pres¬ 
ence of asteroseismic constraints the impact on age estimates of 
a change by +1 in AY/ AZ is indeed negligible, while the contrary 
occurs if mass and radius are available (IValle et al.!l2015l) . 


Because of its potential impact in reducing the number of es¬ 
timated young bulge stars, the effect of a large increase in the 
initial helium content is worthy of a detailed quantitative anal¬ 
ysis. Therefore, we computed a non-standard grid of models 
with the same parameters or input physics and the same values 
of the metallicity Z as in the standard grid, but by significantly 
increasing the helium-to-metal enrichment ratio, i.e. AY/AZ = 5 
instead of 2. Then we repeated the age estimate by adopting this 
helium-enhanced grid in the SCEPtER pipeline. To account for 
possible correlations among the estimates given by the standard 
and helium-enhanced grids, some care was needed in the com¬ 
parison. In fact, the naive differences of the two estimates (and 
their errors) could be obtained by subtracting the estimates pro¬ 
vided by the standard and helium-enhanced grids and by adding 
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Fig. 2. Left : SCEPtER age estimates for the 58 microlensed stars versus the lBensbv et ali d2013l) estimates. Right: difference in age estimates 
between the t wo t echniques versus [Fe/H], A positive value means that the SCEPtER age is greater than the corresponding age estimated by 
iBensbv et al.1 i2013l) . The red dashed line marks the mean difference of -0.08 Gyr between the two sets of age estimates. 




Fig. 4. Left: differences between the age estimates obtained by adopting the grid with AF/AZ = 5 (helium-enhanced) and that with AF/AZ = 2 
(standard). The red dashed line marks the mean of the differences, while the solid red line is a LOESS smoother of the data (see text). Right: 
empirical cumulative distribution functions for the upper limits of the estimated ages obtained by SCEPtER assuming the standard initial helium 
value (black line), the helium-enhanced scenario (blue line), and for the ages estimated bv lBensbv et al.11201 j i) (red line). 


in quadrature the errors in the single estimates. However, since 
such an approach largely overestimates the final errors, we pre¬ 
ferred to follow a more accurate procedure. For each star, we ob¬ 
tained N = 10000 artificially perturbed stars adding a Gaussian 
noise to the observables as detailed in Sect. [2] The artificial star 
ages are then estimated by means of the SCEPtER pipeline with 
both grids. For each of the N artificial stars, we evaluated the 
difference between the age returned by the helium-enhanced and 
the standard grids. The median of these differences was adopted 


as best estimate and the 16th and 84th quantiles were adopted as 
lcr error. 

The results of this analysis are given in Fig. [4] The left panel 
of the figure shows the differences in the estimates computed 
with the helium-enhanced and the standard grid of models. It 
appears that the change in the initial helium content has a lit¬ 
tle influence in age estimates. Adopting the helium-enhanced 
grid, the mean increase in age estimates is 0.35 Gyr (95% con¬ 
fidence interval [0.26; 0.44] Gyr). As expected, the influence of 
the initial helium abundance change is greater at high metallic - 
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ity. The effect is shown in the figure by superimposing a LOESS 
smoothed on the data: the higher the metallicity, the larger the 
discrepancy, until a plateau of about 0.6 Gyr for [Fe/H] > 0. 
However, such an increase in the age estimates is too small to 
rule out the evidence of a young component in the microlensed 
population. 

As one can see in Fig. [2] the error bars on the age estimates 
are relevant. An important point to address is thus whether the 
previous finding might result simply from a fluctuation. Since 
we are mainly interested in checking the robustness of the iden¬ 
tification of a sizeable fraction of young stars, we further anal¬ 
ysed our results by considering as age estimates those provided 
by the upper values of the error bars rather than the mean val¬ 
ues. With such a choice, the probability of obtaining by chance 
an age older than the adopted one is 16%. The right panel 
in Fig. [T] shows the empirical cumulative density function for 
the upper values obtained by SCEPtER assuming standard and 
helium-enhanced grids of models, and bv lBensbv et alj d2013h . 
The three cumulative distributions closely agree up to 12 Gyr. 
All three sets of estimates contain 15 stars younger than 9 Gyr 
(26% of the sample, 95% confidence interval [16%; 39%]). For 
stars younger than 5 Gyr, the estimates assuming standard he¬ 
lium abundance show the presence of seven such objects (12% 
of the sample, 95% confidence interval [5%; 24%]), while there 
are four (7% of the sample, 95% confidence interval [2%; 18%]) 
in the helium-enhanced computations. 

Thus, even in the most conservative scenario in which we 
contemporaneously adopt a high initial helium abundance and 
the upper values of the age estimates, there is strong evidence of 
a young population among the microlensed bulge stars. 


4.2. Effect of neglecting microscopic diffusion 

An important source of systematic bias in the age estimate of 
dwarf stars is the treatment of element microscopic diffusion 
adopted in stellar model computations. Depending on its ef¬ 
ficiency, this process affects the evolution of both the surface 
chemical abundance - and thus [Fe/H] - and the internal struc¬ 
ture of low-mass stars. As a direct consequence, the stellar pa¬ 
rameters determined by comparing theoretical models and obser¬ 
vations depend on the diffusion efficiency adopted. In lValle et alJ 
d2014bid ) we showed that grid-based estimates of mass, radius, 
and age are significantly different between grids of models that 
take into account diffusion and those that neglect it. From the ob¬ 
servational point of view, microscopic diffusion has been proved 
to be efficient in the Sun (see e.g. lBahcall et al.ll2001 ). while its 
actual effi ciency in Galac t ic globular cluster stars i s debated (see 
e.g. Gruvters et al. 2014tlNordlander et alJl2012HGratton et alj 


1201 itlKomelai][2007th 

The lBensbv et alj ( 2013 ) results are based on stellar models 
that consider only the diffusion of helium and not that of metals. 
On the contrary, our standard grid follows the diffusion of helium 
and metals, as mentioned in Sect. 12.11 Since our main aim is to 
verify whether the identification of young stars is robust against 
theoretical biases, we think that it is worth exploring the effects 
of completely neglecting microscopic diffusion in stellar calcu¬ 
lations. Such a choice would in fact maximize age estimates, 
the other parameters being fixed, and consequently minimize 
the sample of young stars. We computed a non-standard grid 


2 A LOESS (LOcal regrESSion) smoother is a non-parametric locally 
weighted polynomial regression technique that is often used to highlight 
the underlying trend of scattered data (see e.g jFeigelson & Bahu|[20 1 2t 
lYenab les & Rio l ev! 20021) . 


of models with exactly the same characteristics as the standard 
model but without the inclusion of a diffusion mechanism. Then, 
we repeated the age estimate by adopting this non-standard grid 
in the SCEPtER pipeline. The differences between the age es¬ 
timates with and without diffusion have been computed in the 
same way as discussed in Sect. l4.H for the helium-rich case. 

The results are shown in the left panel of Fig. [5] As ex¬ 
pected (see e.g. IValle et aPl2014 a). the neglect of diffusion in 
stellar models leads to older ages. However, such an effect gets 
smaller and smaller at decreasing ages owing to the long time 
scale of diffusion processes. The maximum increase of the es¬ 
timated age is about 2.8 Gyr at 8 Gyr and lower than 1.5 Gyr 
for stars younger than 6 Gyr. The mean age increase is indeed 
much smaller being of 0.40 Gyr (95% confidence interval [0.17; 
0.63] Gyr) for stars younger than 5 Gyr. The coloured region 
in the figure highlights stars with estimated ages higher than 
10 Gyr. The decreasing trend in the age differences for ages 
higher than about 10 Gyr is an edge-effect, a common feature of 
any grid -based technique, as extensively discussed in Valle et al.l 
d2014bl lah. Such an effect is the consequence of the lack of mod¬ 
els older than 15 Gyr in the grid, which thus limits the possibility 
of age overestimation for these old stars. 

The effect of the diffusion on age estimates can also be seen 
in the right panel of Fig. [5] which shows the empirical cumula¬ 
tive density functions for the upper values for the age estimates 
obtained considering or neglecting microscopic diffusion. 

It is apparent that even neglecting the microscopic diffusion 
into stellar model computations, there is no evidence against the 
conclusion of the presence of a relevant young stellar population. 
Even in the most unfavourable case in which the non-standard 
grid without diffusion and the upper values of the age estimates 
are used there are 4 (7% of the sample, 95% confidence interval 
[2%; 18%]) younger than 5 Gyr and 11 (19% of the sample, 95% 
confidence interval [10%; 32%]) younger than 9 Gyr. 

4.3. Effect of increasing the superadiabatic convection 
efficiency 

Another source of systematic bias in age determination of dwarf 
stars is the efficiency of convective transport in superadiabatic 
regimes adopted in stellar computations. As is well known, mod¬ 
ern evolutionary codes still have to rely on oversimplified treat¬ 
ments of superadiabatic convection which depend on free param¬ 
eters to be calibrated with observations. The most commonly 
adopted treatment is the mixing-length theory ( Bohm-Vitensel 
119581) . where the efficiency of convective transport depends on 
the free parameter a m i, usually calibrated with the Sun. A solar 
calibrated mixing-length has been adopted both in our standard 
models and in those used by iBensbv et alj (20131 ). However, 
there is no a priori compelling reason to guarantee that a so¬ 
lar calibration also holds in evolutionary phases and/or masses 
different from that of the Sun. Moreover, there are indications 
that a slightly higher value (a m i = 1.90) than the solar cal¬ 
ibrated is in better agreement with the colour-magnitude dia¬ 
grams of old and intermediate age Galactic globular clusters (see 
e.g. lDeirOmodarme et aljr2012l) . For this reason we think that 
it is worth discussing the effect on the age estimates of a change 
in the mixing-length value. As in the two previous tests, we fo¬ 
cused only on a variation that would increase stellar ages, thus 
reducing the number of young bulge stars. 

We computed a non-standard grid of models with the same 
characteristics of the standard one but computed with a m \ = 
1.90. Then, we repeated the age estimate by adopting this non¬ 
standard grid in the SCEPtER pipeline. The differences in ages 
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Fig. 5. Left : Differences in age estimates assuming the grid that neglects microscopic diffusion and the standard one. The solid red line is a 
LOESS smoother of the data (see text). The yellow region shows the stars older than 10 Gyr. Right: empirical cumulative distribution functions 
for the upper limits of the estimated ages obtained by SCEPtER adopting the grid with element diffusion (black line), and the grid that neglects it 
(red line). 


between the a m \ = 1.90 and the ar m i = 1.74 (solar calibrated) 
cases have been computed as in the previous tests (see Sects. |4~TI 
and 14.21) . The results are presented in Fig. [6] The age increase 
due to the higher mixing-length value is lower than the increase 
due to neglecting microscopic diffusion, and its maximum is 
about 1.5 Gyr. For stars younger than 5 Gyr, the mean increase 
in age estimate is 0.54 Gyr (95% confidence interval [0.28; 0.81] 
Gyr). 

The right panel of Fig.[6]shows the empirical cumulative den¬ 
sity functions for the upper values of the age estimates obtained 
with grids of stellar models with the two considered values of 
a m \. The differences between the two distributions are small. 
Even in the most unfavourable case, i.e. simultaneously adopt¬ 
ing the grid with ar m i = 1.90 and the upper values of the age, 
estimates provide 4 (7% of the sample, 95% confidence interval 
[2%; 18%]) younger than 5 Gyr and 13 (22% of the sample, 95% 
confidence interval [13%; 36%]) younger than 9 Gyr. 

In conclusion, even accounting for a higher mixing-length 
value, the presence of a relevant young stellar population is con¬ 
firmed. 


5. Conclusions 

We independently verified the results by iBensbv et alJ ( 2013 ) 
on the age of 58 microlensed dwarf and subgiant stars in the 
Galactic bulge by means of completely different stellar models 
and estimation techniques. 

To this purpose we used the SCEPtER pipeline dValle et al.l 
l2014biah . which allows mass, radius, and age of single stars to be 
estimated by relying on a fine grid of precomputed stellar tracks 
and a set of observational constraints. Here we adopted the ef¬ 
fective temperature the surface gravity, and [Fe/H] provided by 
IBensbv et al.ld2013 ). The standard grid was composed by mod¬ 
els from ZAMS up to the helium flash at the RGB tip, for mass 
in the range [0.60; 1.60] M 0 . 


The performance of the estimation technique for the con¬ 
sidered evolutionary phases was assessed by means of a Monte 
Carlo simulation. The results of this test showed that the adopted 
method is unbiased and it can be safely employed for the re¬ 
quired age estimations. 

We confirmed the results by IBensbv et alJ d2013l) . computed 
adopting a different estimation technique and a different grid of 
stellar models. This confirmation adds considerable robustness 
to the age estimates and in particular to the discovery of a size¬ 
able fraction of young bulge stars. We found in fact that the 
analysed sample of microlensed stars contains 16 stars (about 
one third of the sample) younger than 5 Gyr and 27 (about one 
half of the sample) younger than 9 Gyr. Even conservatively ac¬ 
counting for the errors in the obtained stellar ages, assuming the 
upper values as best age estimates, we found that 15 stars were 
definitely younger than 9 Gyr. 


Given the importance of the result and its widespread impli¬ 
cations, we analysed in detail possible biases in the age estimates 
due to some assumptions adopted in computing the grid of stellar 
models. We first checked if helium rich models would be able to 
rule out the evidence of the presence of young and intermediate 
age stars as suggested bv lNataf & Gouldl ( 2012 ). To this purpose 
we repeated the age estimation procedure relying on a helium- 
enhanced grid of models computed on purpose by adopting a 
helium-to-metal enrichment ratio AY/AZ - 5. The result was a 
median increase in the age of the stars with [Fe/H]> 0 by about 
0.6 Gyr, largely insufficient to reconcile the stellar age estimates 
with the picture of a totally old bulge suggested by photometric 
studies of the main-sequence turnoff (e.g. IZoccali et ak 20031 : 


iBrown et aljl201(t IClarkson et al.ll2()TTL I Valenti et al.ll2013h . To 

exclude the possibility that the results could arise only by a fluc¬ 
tuation, we also verified the robustness of our finding assuming 
the 1 cr upper values of age estimates as best values rather then 
the mean ones. Even with this extreme choice and adopting the 
helium-enhanced grid of stellar models there is evidence of 4 
and 15 stars younger than 5 Gyr and 9 Gyr, respectively. 
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Fig. 6. Left : same as in the left panel of Fig. [5] but for age differences due to assuming the grid with a ml = 1.90 and the reference case with 
Ofmi = 1.74. Right : empirical cumulative distribution functions for the upper limits of the estimated ages obtained by SCEPtER adopting the 
standard grid with or ml = 1.74 (black line), and that with or ml = 1.90 (red line). 


Another important source of systematic bias in the age esti- 
mate of dwarf stars is the treatment of element diffusion (Valle 


et al. 2014biah . We verified that even neglecting the diffusion in 
the stellar models adopted in the estimate procedure is not able to 
increase the age enough to rule out the presence of a young and 
intermediate age population. We obtained a maximum increase 
of about 1 Gyr for stars younger than 5 Gyr and of about 2.8 Gyr 
for stars with ages of about 8 Gyr. The mean age increase is in¬ 
deed much smaller, about 0.40 Gyr for stars younger than 5 Gyr. 
Even in the most unfavourable case in which the upper values of 
the age estimates are used, the number of stars younger than 5 
and 9 Gyr are, respectively, 4 and 11. 

Finally, we checked the effect of changing the calibration 
of the mixing-length from a m \ = 1.74 (our solar value) to 
a m i = 1.90 - a value calibrated on the RGB colour of old and 
intermediate age clusters - in the stellar models used to estimate 
the age. As in the previous cases, the resulting age increase is not 
old enough to reject in a statistically meaningful way the pres¬ 
ence of young bulge stars. In fact, we found a maximum age 
increase of about 1.5 Gyr, and a mean increase of 0.54 Gyr for 
stars younger than 5 Gyr. Even adopting the upper values of the 
age estimates, the stars younger than 5 and 9 Gyr are, respec¬ 
tively, 4 and 13. 

In conclusion, the presence of a young and intermediate age 
population among the microlensed bulge appears very robust and 
statistically confirmed. 
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Appendix A: Expected biases on grid estimates 

We present here an evaluation of the expected errors and possible 
bias of the technique adopted for age estimates. This informa¬ 
tion helps to evaluate the accuracy of the results presented in the 
paper. 

For this test we sampled N - 50 000 artificial stars from the 
recovery grid. To mimic the microlensed sample characteristics, 
we restricted the sampling to stars up to the very first phase of 
RGB. The observational parameters of these stars are then sub¬ 
jected to a random Gaussian perturbation, assuming as observa¬ 
tional errors 100K in T e ff, 0.2 dexin logg, and 0.1 dex in [Fe/H]. 
The star masses and ages were then estimated by means of the 
SCEPtER technique. 

The resulting relative errors in mass and age are presented in 
Fig. IA.ll where a positive value in the relative error implies that 
the technique overestimates the considered parameter. The fig¬ 
ure shows the trend of the relative errors versus the true mass of 
the star and its relative agfl To better follow the resulting trends 
in the figure we superimpose the median and the hr error enve¬ 
lope, computed over a moving window (see I Valle et alJ[2014aL 
for details on the envelope computation). The figure shows some 
clear signatures of edge effects, such as the trend in the rela¬ 
tive error of mass estimate versus the true mass (top left panel 
in Fig. IA. II) or the huge distortion of age estimate at low rela¬ 
tive age (bottom right panel in the figure), extensively discussed 
in IValle et all d2014blah . As an example, the trend in the top 
left panel of the figure at the high mass edge (M = 1.6 M Q ) 
originates from the impossibility of obtaining a mass estimate 
above 1.6 M e , the maximum value considered in the grid con¬ 
struction. Similar considerations apply to the other grid edges 
(IValle et al.ll2014al) . The tail toward age overestimation, visible 
at low masses in the bottom left panel, is due to the fact that for 
stars of mass lower than about 0.7 M 0 the grid does not contain 
models at high relative age because of the constraint of the 15 
Gyr maximum age mentioned in Sect. 12.11 In other words for 
these masses the grid does not reach the later central hydrogen 
burning phases, which leads to more precise age estimates (see 
bottom right panel in Fig. IA.ll ). The lack of these models causes 
the asymmetry in the envelope, typical of the first evolutionary 
phases. As expected (see e.g. lGai et all201 il) . mass and age es¬ 
timates become more difficult in the subgiant phase (relative age 
higher than 1), since tracks are more closely packed in this grid 
zone. Moreover, for relative ages higher than about 1.1a slight 
mass underestimation distortion appears, with a corresponding 
increase in age estimate (see the lower envelope boundary in the 

3 The relative age is defined as the ratio between the current age of the 
star and its age at the central hydrogen depletion. It is conventionally 
set to 0 at the ZAMS. 


top right panel of Fig. [AT] and the corresponding upper bound¬ 
ary in the bottom right panel). The trend is due to the packing 
of the tracks toward the RGB ascension; low mass tracks pack 
more strictly than those of high mass, therefore the mass is more 
easily underestimated than overestimated. 

For our purposes, the most interesting feature is the appar¬ 
ent lack of bias in the median of mass and age estimates visible 
in Fig. IA. ll fwith the exception of the first 20% of the evolution 
where edge effects dominate). It also appears from the bottom 
right panel of the figure that in the subgiant phase the age esti¬ 
mates present a tail toward overestimation. 

It is often stated in the literature that mass and age estimates 
computed neglecting the evolutionary time step suffer from a 
bias and they could be improved by account i ng for this correc¬ 
tion ( see e.g. iJdreenserT & Lindegren 2005; Casagr ande et all 
1201 ill . This statement also appears in Bensbv et alJ ( 2013 ) in 
relation to the microlensed star age estimates. However, this 
is not always true since the effect is strictly related to the em¬ 
ployed technique and the considered evolutionary stages. In par¬ 
ticular, for the results presented in this paper, the weighted es¬ 
timates obtained adopting the evolutionary time step as weight 
are clearly biased toward age overestimation (bottom right panel 
in Fig. IA, I b . The effect is most evident in the subgiant phase, 
where the bias is higher than 100%. In this phase, where tracks 
are closely crowded, the time step weight can make a large dif¬ 
ference, leading to the selection of low-mass high-age models. 
Figure IaTI shows these trends by reporting the lcr envelope of 
the weighted estimates. Therefore, in the light of this result, in 
the present work we adopt unweighted estimates as a reference 


scenario. 
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Fig. A.l. Top row : left Monte Carlo relative errors on mass estimate as a function of the true mass of the stars. The black lines show the position 
of the median (dashed line) and the lcr envelope (solid line). The red lines show the corresponding quantities, but for weighted estimates. Right : 
same as in the left panel , but as a function of the relative age of the stars. Bottom row : same as the top row , but for age estimates. 
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Table 1. Ages and masses determined by the SCEPtER pipeline. 


Object 

age 

(Gyr) 

—lcr 
(Gyr) 

lcr 

(Gyr) 

mass 

(Mo) 

-lcr 

(M 0 ) 

lcr 

(M Q ) 

[Fe/H] 

A [Fe/H] 

OGLE-2012-BLG-1156 

11.1 

6.5 

3.2 

0.76 

0.05 

0.09 

-1.89 

0.25 

OGLE-201 l-BLG-0969 

10.4 

5.0 

3.8 

0.79 

0.05 

0.13 

-1.57 

0.25 

MOA-2010-BLG-285 

12.0 

3.9 

2.5 

0.79 

0.04 

0.07 

-1.21 

0.10 

MOA-2010-BLG-07 8 

13.7 

8.4 

0.8 

0.82 

0.04 

0.25 

-1.00 

0.34 

MOA-2011-BLG-104 

10.2 

5.4 

3.9 

0.83 

0.05 

0.13 

-0.85 

0.15 

OGLE-2012-BLG-0270 

10.8 

7.0 

3.4 

0.81 

0.04 

0.07 

-0.84 

0.13 

MOA-2012-BLG-187 

12.7 

3.1 

1.9 

0.81 

0.03 

0.05 

-0.74 

0.08 

MOA-2009-BLG-493 

13.7 

9.4 

1.0 

0.72 

0.04 

0.08 

-0.74 

0.15 

OGLE-2009-BLG-076 

13.0 

4.9 

1.6 

0.80 

0.04 

0.04 

-0.72 

0.11 

MOA-2009-BLG-133 

14.5 

6.9 

0.3 

0.74 

0.02 

0.06 

-0.69 

0.07 

OGLE-2012-BLG-0563 

10.7 

6.4 

3.2 

0.81 

0.03 

0.04 

-0.66 

0.07 

OGLE-2012-BLG-1279 

13.4 

5.1 

1.4 

0.78 

0.02 

0.03 

-0.62 

0.06 

MOA-2010-BLG-167 

14.1 

1.6 

0.1 

0.82 

0.00 

0.03 

-0.60 

0.05 

MOA-2012-BLG-532 

8.9 

6.6 

5.0 

0.90 

0.10 

0.45 

-0.55 

0.21 

MOA-2009-BLG-475 

7.9 

6.7 

5.8 

0.84 

0.06 

0.08 

-0.52 

0.20 

MACH-1999-BLG-022 

14.1 

3.9 

0.7 

0.83 

0.04 

0.08 

-0.49 

0.17 

OGLE-2012-BLG-1217 

12.2 

3.4 

2.2 

0.83 

0.03 

0.04 

-0.41 

0.07 

MOA-2010-BLG-049 

14.0 

2.8 

0.8 

0.84 

0.02 

0.06 

-0.40 

0.07 

MOA-2010-BLG-446 

0.9 

0.6 

2.6 

1.01 

0.05 

0.05 

-0.40 

0.08 

OGLE-2008-BLG-209 

11.3 

5.4 

3.0 

0.90 

0.06 

0.18 

-0.31 

0.09 

MOA-201 l-BLG-090 

14.5 

0.5 

0.4 

0.86 

0.02 

0.00 

-0.26 

0.05 

OGLE-2012-BLG-1526 

13.2 

5.1 

1.4 

0.88 

0.02 

0.13 

-0.24 

0.06 

MOA-2012-BLG-391 

10.0 

4.2 

4.1 

0.95 

0.09 

0.15 

-0.24 

0.07 

MOA-2009-BLG-489 

14.3 

2.5 

0.6 

0.86 

0.02 

0.05 

-0.21 

0.07 

MOA-2011-BLG-174 

2.5 

2.1 

2.5 

1.05 

0.05 

0.07 

-0.18 

0.09 

OGLE-2012-BLG-1534 

7.3 

1.0 

1.1 

1.04 

0.05 

0.05 

-0.15 

0.04 

MOA-2012-BLG-202 

14.6 

0.4 

0.2 

0.88 

0.02 

0.02 

-0.15 

0.13 

OGLE-2012-BLG-0617 

14.2 

4.5 

0.5 

0.89 

0.03 

0.10 

-0.14 

0.09 

MOA-2012-BLG-410 

8.9 

3.6 

4.2 

0.99 

0.10 

0.15 

-0.14 

0.07 

OGLE-2012-BLG-0816 

1.8 

0.2 

2.4 

1.57 

0.34 

0.03 

-0.10 

0.14 

OGLE-2012-BLG-0211 

9.9 

2.9 

2.9 

0.99 

0.08 

0.10 

-0.03 

0.08 

MOA-201 l-BLG-234 

6.2 

4.9 

3.2 

0.96 

0.04 

0.05 

-0.02 

0.08 

OGLE-2011-BLG-1105 

3.1 

2.6 

6.9 

0.93 

0.06 

0.07 

0.00 

0.17 

MOA-2010-BLG-523 

2.8 

1.0 

4.2 

1.42 

0.33 

0.17 

0.06 

0.14 

OGLE-2012-BLG-1274 

8.5 

1.0 

1.1 

1.03 

0.04 

0.04 

0.07 

0.04 

OGLE-2012-BLG-0521 

6.8 

4.0 

6.3 

1.12 

0.19 

0.33 

0.09 

0.15 

MOA-201 l-BLG-034 

12.4 

3.2 

2.2 

0.94 

0.05 

0.09 

0.11 

0.15 

MOA-2009-BLG-174 

7.8 

5.0 

3.1 

0.94 

0.04 

0.04 

0.11 

0.08 

MOA-2009-BLG-456 

9.0 

1.9 

2.0 

1.00 

0.06 

0.06 

0.13 

0.10 

MOA-2012-BLG-291 

14.5 

2.8 

0.3 

0.92 

0.06 

0.06 

0.16 

0.24 

OGLE-201 l-BLG-1410 

3.9 

1.7 

9.5 

1.35 

0.41 

0.23 

0.22 

0.37 

MOA-2011-BLG-191 

4.6 

2.3 

3.5 

1.27 

0.18 

0.26 

0.26 

0.14 

MOA-201 l-BLG-445 

4.8 

2.6 

8.5 

1.26 

0.32 

0.31 

0.26 

0.31 

OGLE-2007-BLG-514 

6.6 

2.7 

3.1 

1.09 

0.11 

0.16 

0.29 

0.23 

OGLE-201 l-BLG-0950 

2.2 

1.3 

1.0 

1.29 

0.09 

0.12 

0.33 

0.10 

MOA-2009-BLG-259 

2.2 

0.2 

1.7 

1.59 

0.24 

0.01 

0.34 

0.19 

MOA-2008-BLG-311 

0.3 

0.1 

1.1 

1.14 

0.04 

0.04 

0.35 

0.08 

OGLE-2011-BLG-1072 

5.7 

2.2 

2.5 

1.20 

0.12 

0.17 

0.36 

0.12 

MOA-201 l-BLG-058 

11.0 

4.9 

3.3 

1.00 

0.08 

0.18 

0.37 

0.25 

MOA-2008-BLG-310 

6.2 

1.7 

1.8 

1.09 

0.06 

0.07 

0.41 

0.11 

MOA-2012-BLG-022 

3.5 

2.4 

1.7 

1.15 

0.06 

0.07 

0.42 

0.09 

OGLE-2007-BLG-349 

13.8 

4.1 

0.9 

0.93 

0.05 

0.06 

0.42 

0.26 

MOA-2006-BLG-099 

0.6 

0.3 

3.6 

1.08 

0.07 

0.06 

0.44 

0.21 

OGLE-2006-BLG-265 

9.4 

2.6 

2.8 

1.00 

0.06 

0.07 

0.46 

0.19 

OGLE-2012-BLG-0026 

3.8 

1.6 

9.2 

1.37 

0.41 

0.21 

0.50 

0.44 

MOA-2010-BLG-311 

3.9 

1.7 

2.9 

1.33 

0.19 

0.24 

0.51 

0.19 

MOA-201 l-BLG-278 

8.4 

5.4 

5.2 

1.04 

0.12 

0.33 

0.52 

0.39 

MOA-2010-BLG-037 

3.8 

1.6 

1.7 

1.31 

0.14 

0.23 

0.55 

0.20 


Notes. Column 1: star identifier; Col. 2: stellar age; Cols. 3 and 4: lcr lower and upper uncertainties for age estimates; Col. 5: stellar mass; Cols. 
6 and 7: lcr lower and upper uncertainties for mass estimates; Cols. 8 and 9: [Fe/H] and its error from iBensbv et al J ( 2013 ). 
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